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Abstract 


We  consider  the  approximation  algorithm  of  Leighton  et.  al.  [8]  for  the  multicom¬ 
modity  flow  problem.  We  give  a  more  natural  randomization  strategy  that  is  simpler 
than  the  one  in  [8]  and  results  in  a  better  running  time.  This  strategy  also  applies  to 
several  related  algorithms. 


1  Introduction 

The  multicommodity  flow  problem  is  as  follows:  given  a  network  with  capacity  constraints 
and  commodity  supplies  and  demands,  find  a  flow  that  satisfies  the  demands  without  vio¬ 
lating  capacity  constraints.  The  multicommodity  flow  problem  is  a  classical  problem  that 
has  numerous  applications.  The  concurrent  flow  problem  is  an  optimization  version  of  the 
multicommodity  flow  problem,  where  the  goal  is  to  maximize  the  fraction  of  the  satisfied 
demands,  i.e.,  to  maximize  z  such  that  if  the  demands  are  multiplied  by  z ,  the  resulting 
multicommodity  flow  problem  is  feasible. 

In  this  paper  we  work  with  the  concurrent  flow  problem.  We  denote  the  number  of 
nodes  in  the  network  by  n,  the  number  of  arcs  by  m,  the  number  of  commodities  by  k,  the 
largest  capacity  by  U,  the  largest  demand  by  D,  and  assume  that  capacities  and  demands 
are  integral. 

The  only  known  polynomial  time  algorithms  for  the  problem  are  based  on  polynomial¬ 
time  methods  for  linear  programming,  either  the  ellipsoid  method  [6]  or  the  interior-point 
method  [5].  The  fastest  currently  known  algorithm,  due  to  Vaidya  [12],  takes  advantage 
of  the  matrix  structure  to  achieve  improvement  over  the  underlying  interior- point  method. 
This  algorithm  runs  in  0(fc3-5n3m°-5  log (nDU))  time. 

For  many  applications  it  is  sufficient  to  find  an  approximate  solution  to  the  problem,  i.e., 
a  feasible  solution  that  ships  (l-e)z*  fraction  of  the  demands,  where  z *  is  the  optimal  value. 
In  these  applications,  €  is  a  constant  or  a  slowly  decreasing  function  of  n  (e.g.  £).  Shahrokhi 
and  Matula  [11]  gave  an  approximation  algorithm  for  a  special  case  of  the  concurrent  flow 
problem  with  uniform  capacities.  Their  work  motivated  several  other  papers.  A  faster 
algorithm  for  the  uniform  capacity  case  was  given  by  Klein  et.  al.  [7],  An  algorithm  for 
the  general  version  of  the  concurrent  flow  problem  was  given  by  Leighton  at.  al.  [8].  The 
deterministic  version  of  this  algorithm  runs  in  0(k2(e~ 2  log  j  +  lognlogfc))  minimum-cost 
flow  computations  and  the  randomized  version  in  0(k(e~3  log  j  -flog  n  log  k))  minimum-cost 
flow  computations.  (For  a  survey  of  minimum-cost  flow  algorithms,  see  e.g.  [2].) 
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In  this  paper  we  give  a  modification  of  the  relaxed  optimality  conditions  that  leads 
to  a  slightly  simpler  analysis  of  the  algorithm  and  allows  us  to  show  that  the  simplest 
randomization  strategy,  which  selects  a  commodity  to  be  updated  during  the  next  iteration 
uniformly  at  random,  works  better  than  the  more  complicated  strategy  used  in  [8).  We  show 
an  0(k(e~2  log  £  +  logo  log  A:))  bound  on  the  number  of  minimum-cost  flow  computations. 

From  the  theoretical  viewpoint,  the  dependence  on  e  is  a  major  drawback  of  the  algo¬ 
rithm  of  [8]  since  exponential  precision  is  required  to  solve  the  multicommodity  flow  problem 
exactly.  We  reduce  this  dependence  by  a  factor  of  e-1.  Our  randomization  strategy  appears 
to  be  a  better  choice  in  practice  as  well. 

Our  results  can  also  be  applied  in  a  straightforward  way  to  improve  the  results  of 
[7]  on  the  uniform  multicommodity  flow  problem,  as  well  as  in  the  generalizations  of  the 
multicommodity  flow  algorithm  to  other  linear  programming  problems  [10]. 

A  recent  paper  of  Grigoriadis  and  Khachiyan  [4]  extends  the  results  of  [8]  to  block- 
structured  convex  linear  programs.  This  paper  uses  a  randomization  strategy  similar  to 
ours  in  a  somewhat  different  framework. 

This  paper  is  organized  as  follows.  Section  2  gives  definitions  and  notation  used  in  the 
paper.  Section  3  describes  relaxed  optimality  conditions.  Section  4  discusses  the  exponen¬ 
tial  length  function  and  its  properties.  Section  5  describes  the  algorithm  and  Section  6 
analyses  it.  Section  7  discusses  some  aspects  of  the  algorithm  and  general  applicability  of 
our  techniques. 


2  Definitions  and  Notation 

In  this  paper  we  consider  a  directed  version  of  the  multicommodity  flow  problem;  however, 
the  results  also  hold  for  the  undirected  version  with  straightforward  modifications.  An  input 
to  the  multicommodity  flow  problem  is  a  graph  G  =  (V,  E),  a  capacity  function  u  :  V  — > 
R+,1  and  a  demand  specification  for  commodities.  For  each  commodity  *:!<*<  k ,  the 
specification  contains  the  source  s,,  the  sink  f;,  and  a  nonnegative  demand  d{.  Without 
loss  of  generality  we  assume  that  m  >  n  and  that  G  has  no  parallel  arcs.  We  assume  that 
capacities  and  demands  are  reals  and  denote  the  biggest  capacity  by  U  and  the  biggest 
demand  by  D. 

1 R*  denotes  the  set  of  non  negative  reals. 
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A  multicommodity  flow  /  is  given  by  a  set  of  functions  fl, 

Each  function  fl  must  satisfy  conservation  constraints 

di  if  v  = 

Vv  €  V,  £  fi(v,w)=<  -di  \{v  =  si, 

(u,t/)g  £  (v,tu)gE  0  otherwise. 

We  define  f(v,w)  =  Za<«<fe  fi{v,w).  A  multicommodity  flow  is  feasible  if  the  capacity 
constraints  are  satisfied: 

V(v,  w)  €  £,  /(«,  u>)  <  u(v,  u>). 

The  concurrent  flow  problem  [11]  is  an  optimization  version  of  the  multicommodity  flow 
problem  where  the  objective  is  to  maximize  z  such  that  the  problem  with  demands  zd,  is 
feasible.  An  equivalent  problem  is  to  minimize  A  such  that  the  problem  with  demands  d, 
and  capacities  Au  is  feasible.  Let  A*  denote  the  optimal  value  of  A.  We  define  the  congestion 
by  A(v,  w)  =  f(v,tv)/u(v,w),  and  let  A  =  max^jgg  X(v,w).  A  multicommodity  flow  /  is 
(-optimal  if  A  <  (1  +  t)A*.  In  this  paper  we  consider  the  problem  of  finding  an  e-optimal 
solution  to  the  problem;  throughout  the  paper  we  assume  that  0  <  t  <  l. 

Next  we  introduce  the  length  function,  which  has  its  roots  in  linear  programming  duality. 
The  length  function  l :  E  -*  R+  is  a  nonzero,  nonnegative  function.  Let  dist((v,  w)  denote 
the  distance  in  G  from  v  to  w  with  respect  to  l. 

To  simplify  notation,  we  sometimes  view  length,  capacity,  and  flow  functions  as  vectors 
indexed  by  arcs.  Given  two  vectors  a  and  6,  let  a  o  b  =  Yl(v,w)eE  a(u>  w)b{v,  w). 

Theorem  2.1  For  a  multicommodity  flow  /  and  a  length  function  l  we  have 

k 

A t(v,  w)  o  u(v,  w)  >  distf(si,  t,)d,. 

i=i 

Given  a  length  function  l ,  we  define  the  cost  of  the  flow  of  commodity  i  by  Ci  =  /,  o  l. 
Let  C*(A)  be  the  value  of  the  minimum-cost  flow  fl  satisfying  the  demands  of  commodity 
i  with  costs  l  and  capacities  Au. 

Lemma  2.2  [8]  For  /,  C,-,  and  C*( A)  as  above, 

«oU>£c,>£c-(A). 

i=l  i=l 
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A  multicommodity  flow  /  minimizes  A  iff  there  is  a  length  function  t  for  which  the  above  terms 
are  equal. 


Lemma  2.3  [8]  For  /,  l  as  above, 

EL  QB 

uol 

is  a  lower  bound  on  A*. 

3  Relaxed  Optimality 

The  relaxed  optimality  conditions  are  as  follows: 

V(v,  w)  €  E  (1  +  e)f(v,w)  >  Au(u,  w)  or  u(v,w)l(v,  w)  <  —  (tiof).  (1) 

m 

(l-2f)£C,<£C;(A).  (2) 

i  I 

The  first  condition  is  the  same  as  the  one  in  [8];  it  states  that  either  an  arc  is  close  to  being 
saturated,  or  its  “volume”  is  small  compared  to  the  total  “volume”.  The  second  condition 
states  that  the  total  cost  of  /  with  respect  to  l  is  close  to  the  optimal. 

Theorem  3.1  If  /,  l,  and  c  satisfy  the  relaxed  optimality  conditions,  then  A  is  at  most  (1  + 

0(c))A*. 

Proof.  From  condition  (1),  we  have 

T:  ((1  +  ( )/{v ,  w)i(v,w)  +  A— u  o  l)  >  Xu  o  l. 

(v,w)€E 

Rewriting,  we  get 

(1  +  f)/  o  t  >  Au  o  f(l  -  e) 
or 

i^A«°f</o<  =  ]£Ci<r;i^£ClW  (3) 
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By  Lemma  2.3  and  the  previous  inequality, 


A*  > 


ZiC;(\) 

U  0  / 


(1-0(1-20 
1  +  e 


A, 


and  the  theorem  follows.  ■ 


4  Exponential  Length  Function 


Shahrokhi  and  Matula  [11]  were  the  first  to  use  the  length  function  that  exponentially 
depends  on  flow  in  the  context  of  the  uniform  capacities.  Leighton  et.  ai.  introduced  the 
following  generalization  of  this  idea. 

Define 


l(v,w) 


ga\(v,tv) 

u(v,w) 


(4) 


Lemma  4.1  [8]  Suppose  /  is  defined  by  (4)  and  a  >  (1  +  0A  1  ln(77ie-1).  Then  the  first 
relaxed  optimality  condition  (1)  is  satisfied. 

Proof.  Suppose  (1  4-  e)f(v,w)  <  Xu(v,w).  Then  A(v,  w)  <  A/(l  +  c).  We  need  to  show 
that  u(v,  w)((v,  w)  <  o  l. 

Note  that 

-  uol  =  —Y  eQ^v'w)  >  —eaX. 
m  m  e  m 

To  obtain  the  last  inequality,  we  bounded  the  sum  of  nonnegative  terms  by  the  largest  one. 
On  the  other  hand, 

u(v,w)l(v,w)  =  eoA<v'u,)  <  e°T77  =  eoAe^. 

By  the  assumption  of  the  lemma,  the  second  term  is  at  most  e- and  the  lamma 
follows.  ■ 


5  Algorithm  Description 

The  algorithm  maintains  a  multicommodity  flow  /  and  the  corresponding  exponential  length 
function  l  such  that  /  satisfies  the  demands.  By  Lemma  4.1,  the  first  relaxed  optimality 
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condition  is  always  satisfied.  While  the  second  condition  is  not  satisfied,  the  algorithm  picks 
a  random  commodity  and,  if  the  potential  function  uol  decreases,  replaces  a  fraction  of  its 
flow  by  the  minimum-cost  flow  with  respect  to  the  cost  function  l.  The  expected  decrease 
in  the  potential  function  is  large,  ensuring  quick  termination  of  the  algorithm  with  high 
probability. 

Finding  initial  solution  The  algorithm  starts  with  an  initial  solution  /  such  that  A  < 
k\*.  Such  a  solution  is  obtained  by  finding,  for  each  commodity  i,  a  maximum  flow  <7,  in  the 
network  with  source  Si,  sink  t;,  and  capacity  function  u,  and  setting  fi(v,w)  =  <7i(u,uj)j^|, 

where  |<ft|  denotes  the  value  of  <7,.  It  is  easy  to  see  that  for  the  resulting  flow  /,  A  <  kXm. 
The  length  function  l  defined  by  the  initial  flow  is  also  computed. 

Improving  the  current  solution  At  each  iteration,  the  algorithm  iteratively  improves 
the  current  flow,  until  the  desired  approximation  precision  is  obtained.  An  iteration  ( rerout¬ 
ing  step)  works  as  follows. 

1.  Select  commodity  i  from  the  set  {1,...,/;}  uniformly  at  random. 

2.  Consider  the  network  with  capacity  function  Au  and  cost  function  t.  Compute  the 
minimum  cost  flow  /*  in  this  network  satisfying  demands  for  the  commodity  t. 

3.  Define  f-  =  (1  —  a)fx  +  <7/*,  and  let  f  be  the  multicommodity  flow  obtained  from  / 
by  replacing  /;  by  /'. 

4.  Compute  the  length  function  l'  for  /'. 

5.  Ifuo/>uo£'  then  replace  /  by  /'. 

The  parameters  a  and  a  are  set  as  follows:  a  —  2(  1  +  c)A-1e-1  ln(me-1)  and  a  -  The 
values  of  these  parameters  are  not  updated  at  every  iteration,  but  the  values  are  updated 
when  A  decreases  by  at  least  a  factor  of  two  from  its  value  during  the  time  of  the  last  update. 
Intuitively,  an  iteration  of  the  algorithm  replaces  a  o  fraction  of  a  random  commodity  by 
the  same  fraction  of  the  optimal  flow  of  this  commodity.  The  flow  is  updated  only  if  the 
potential  function  no  i  decreases. 


Nonscaling  and  scaling  algorithms  A  simple  implementation  of  the  algorithm  is  to 
set  e  to  the  desired  value  at  the  very  beginning,  and  improve  the  flow  until  it  becomes 
c-optimal. 

An  alternative  is  to  scale  e.  Note,  however,  that  we  assumed  that  e  <  1.  Thus  the  scaling 
algorithm  starts  with  a  constant  <  (say  1)  and  finds  an  e-optimal  solution  /(0).  Then  at 
jth  scaling  iteration,  e  is  reduced  by  a  factor  of  two  and  a  new  e-optimal  solution  /<■*)  is 
computed  starting  from  until  the  desired  precision  is  reached. 

As  we  shall  see  later,  the  use  of  scaling  improves  the  running  time  of  the  algorithm  by 
a  factor  of  e_1  if  e  is  small  enough. 

Termination  detection  The  above  description  assumes  that  we  know  when  the  current 
flow  becomes  e-optimal.  This  is  not  the  case;  however,  we  can  test  for  e-optimality  by 
computing  minimum-cost  flows  of  all  commodities.  This  can  be  done  every  k  iterations 
while  increasing  the  worst-case  running  time  bound  by  a  constant  factor. 


6  Algorithm  Analysis 


For  the  purpose  of  the  analysis,  we  define  the  potential  function  $  by  $  =  u  o  i . 

The  next  lemma  shows  that  when  fi  is  rerouted,  $  decreases  by  almost  aa(C,  -  C’( A)). 
The  proof  of  this  lemma  is  similar  to  the  proof  of  [8]  showing  that  $  decreases  significantly 
if  a  “bad”  commodity  is  rerouted. 

Let  $  and  <&„  be  the  values  of  the  potential  function  before  and  after  the  rerouting, 
respectively. 


Lemma  6.1 

*  -  $„  >  a<7(C,  -  C;(A)  -  €C<). 

Proof.  Let  l  and  ia  be  the  length  functions  before  and  after  the  rerouting.  By  Taylor’s 
theorem,  for  |f|  <  c/4  <  1/4  we  have 

ex+t  <  er  +  tex  +  ^|t|ex. 


Thus 


£a(v,w)  <  l(v,w)  + 


aa 


u(v,w) 


(/*(u,  w)  -  fi(v,  w))t(v,  w)  + 


CQ<7 


2u(v,tn) 


\f-(v,w)- fi(v,w)\l(v,w) 
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(note  that  f*(v,w )  -  /,(v,  u>)  <  A u(v,w),  and  recall  the  choice  of  a).  We  have 


t 


$  -  $o  =  {t-  la)oU 

>  aa(Ci  -  C*(A))  —  eacrCi. 

The  last  line  follows  by  definition  of  C{  and  C*(A)  and  c y  the  fact  that  C{  >  C*( A)  >  0.  ■ 


Note  that  if  the  current  flow  does  not  satisfy  the  second  relaxed  optimality  condition, 
then  the  expected  value  of  (7,  —  C*( A)  is  large  and  $  decreases  significantly.  The  following 
lemma  formalizes  this  statement. 

Lemma  6.2  Suppose  <  a  <  ^  and  /  does  not  satisfy  the  second  relaxed  optimality 

2 

condition.  Then  the  expected  decrease  in  $  due  to  a  rerouting  step  is 


Proof.  Since  the  algorithm  selects  a  commodity  uniformly  at  random,  we  have 


£[*-*.]  > 


> 


£a<r£(C,- -  c;( A)-«Ci) 

i 

sjk'E«i-K)Ci-cn+'‘Eci) 


aa  1  -  e 

Ti  +  t 

«(T). 


cA$ 


t 


The  third  line  follows  from  the  assumption  that  /  does  not  satisfy  the  second  relaxed 
optimality  condition  and  since  by  (3) 


1  -  f 

1  +  e 


<£c.. 

i 


The  last  line  follows  from  the  assumption  on  a.  ■ 


Note  that  since  a  and  a  are  updated  every  time  A  decreases  by  a  factor  of  two,  conditions 
of  the  lemmas  4.1  and  6.2  are  always  satisfied. 

Next  we  analyze  the  nonscaling  algorithm. 
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Lemma  6.3  Suppose  the  parameters  a  and  a  are  set  when  the  congestion  is  A0.  In  expected 
0(e~3k  log  j)  iterations,  either  the  algorithm  terminates  or  A  decreases  to  Ao/2  or  less. 

Proof.  Initially,  $  <  me4''0  since  the  congestion  on  any  arc  is  at  most  Ao-  Note  that  until 
A  decreases  by  a  factor  of  two,  the  congestion  on  some  arc  is  at  least  Ao/2,  so  $  >  eaA° /2. 
Thus  $  cannot  decrease  by  more  than  a  factor  of  meoA°/2  without  A  getting  below  Ao/2. 

By  Lemma  6.2,  0{c~2k)  iterations  reduce  $  by  a  factor  of  two.  $  can  be  halved  at 
most  0(log(meoAo/2))  =  0(c_1  log  j)  times  before  A  is  halved.  Therefore  in  O(c-3fclog  j) 
iterations  (expected),  the  algorithm  terminates  or  A  is  halved.  ■ 

Note  that  the  work  done  by  an  iteration  of  the  algorithm  is  dominated  by  a  minimum- 
cost  flow  computation.  Combined  with  the  above  lemma  and  the  fact  that  initially  A  <  k A*, 
this  yields  the  following  result. 

Theorem  6.4  The  nonscaling  algorithm  runs  in  expected  O(e~3klog  A: log  ^))  minimum-cost 
flow  computations. 

To  analyze  the  scaling  algorithm,  we  need  the  following  version  of  Lemma  6.3. 

Lemma  6.5  Suppose  the  parameters  a  and  a  are  set  when  the  congestion  is  A0,  and  the  initial 
flow  is  0(c)-optimal.  Then  in  expected  0(e~2k  log  ~)  iterations,  either  the  algorithm  terminates 
or  A  decreases  by  at  least  a  factor  of  two. 

Proof.  Recall  that  e  <  1,  so  Ao  =  0(  A").  The  proof  is  similar  to  that  of  Lemma  6.3,  except 
that  the  upper  bound  on  $  is  meal1+0^^A’^  and  the  lower  bound  is  eaA’,  so  4>  needs  to 
decrease  by  a  factor  of  mea°WXo.  ■ 

The  scaling  algorithm  starts  with  e  =  0(1),  and  obtains  an  O(l)-optimal  solution  in 
O(fclogfclogm)  iterations  by  Theorem  6.4.  Then  the  scaling  process  starts.  Since  the 
number  of  rerouting  steps  needed  to  reduce  e  by  a  factor  of  two  is  proportional  to  €~2,  the 
last  scaling  iteration  dominates.  This  iteration  terminates  in  0(c~2fclog  j)  rerouting  steps. 
We  have 

Theorem  6.6  The  scaling  algorithm  runs  in  expected  0(k(e~ 2  log  j  +  log  k  log  n))  minimum- 
cost  flow  computations. 
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7  Remarks 


First,  we  would  like  to  address  the  selection  of  a  minimum-cost  flow  subroutine.  Since  the 
costs  in  the  minimum-cost  subproblems  can  be  big  due  to  the  use  of  the  exponential  length 
function,  Orlin’s  strongly  polynomial  algorithm  [9]  appears  to  be  the  best  choice.  However, 
it  is  enough  to  solve  the  subproblems  approximately,  so  the  costs  can  be  rounded  to  small 
integers.  When  this  is  done,  the  cost  scaling  algorithm  of  Goldberg  and  Tarjan  [3]  or  the 
double  scaling  algorithm  of  Ahuja  et.  al.  [1]  become  a  better  choice.  The  selection  of  the 
subroutine  is  discussed  in  more  detail  in  [8]. 

The  deterministic  version  of  the  algorithm  of  Leighton  et.  al.  [8]  finds  a  commodity 
i  with  the  biggest  C{  -  C*( A)  and  reroutes  this  commodity.  Our  variant  of  the  relaxed 
optimality  conditions  also  can  be  used  to  analyze  this  algorithm.  Since  the  maximum  of 
a  set  of  numbers  is  at  least  as  big  as  the  average,  the  deterministic  choice  of  i  gives  an 
improvement  in  $  that  is  at  least  as  big  as  the  expected  improvement.  Thus  the  bound 
on  the  number  of  rerouting  steps  apply  for  the  deterministic  version  of  the  algorithm. 
Each  deterministic  rerouting  step,  however,  requires  k  minimum-cost  flow  computations  (to 
compute  C*(A)  for  every  commodity  j)\  in  contrast,  a  randomized  rerouting  step  requires 
a  single  minimum-cost  flow  computation. 

Our  randomization  strategy  can  also  be  applied  to  the  algorithm  of  Klein  et.  al.  [7]. 
This  improves  the  running  time  bound  of  the  randomized  algorithm  of  [7]  by  a  factor  of  <_l . 
We  omit  the  details,  which  are  straightforward  given  those  of  [7],  [8],  and  this  paper.  A 
recent  paper  of  Plotkin,  Shmoys,  and  Tardos  [10],  that  extends  the  results  of  [8]  to  a  more 
general  class  of  linear  programming  problems,  also  extends  our  randomization  strategy  to 
the  more  general  framework. 


Acknowledgements 


The  autl  or  would  line  to  thank  £v a  Tardos  and  Serge  Plotkin  for  helpful  discussions,  and 
Robert  Kennedy  and  Tomasz  Radzik  for  comments  on  a  draft  of  this  paper. 

References 

[1]  R.  K.  Ahuja,  A.  V.  Goldberg,  J.  B.  Orlin,  and  R.  E.  Tarjan.  Finding  Minimum- 
Cost  Flows  by  Double  Scaling.  Technical  Report  STAN-CS-88-1227,  Department  of 


10 


Computer  Science,  Stanford  University,  1988. 


[2]  A.  V.  Goldberg,  E.  Tardos,  and  R.  E.  Tarjan.  Network  Flow  Algorithms.  In  B.  Korte, 
L.  Lovasz,  H.  J.  Promel,  and  A.  Schrijver,  editors,  Flows,  Paths,  and  VLSI  Layout, 
pages  101-164.  Springer  Verlag,  1990. 

[3]  A.  V.  Goldberg  and  R.  E.  Tarjan.  Finding  Minimum-Cost  Circulations  by  Successive 
Approximation.  Math,  of  Oper.  Res.,  15:430-466, 1990.  A  preliminary  version  appeared 
in  Proc.  19th  ACM  Symp.  on  Theory  of  Comp.,  7-18,  1987. 

[4]  M.  D.  Grigoriadis  and  L.  G.  Khachiyan.  Fast  Approximation  Schemes  for  Convex 
Programs  with  Many  Blocks  and  Coupling  Constraints.  Technical  Report  DCS-TR- 
273,  Department  of  Computer  Science,  Rutgers  University,  1991. 

[5]  N.  Karmarkar.  A  New  Polynomial-Time  Algorithm  for  Linear  Programming.  Combi- 
natorica,  4:373-395,  1984. 

[6]  L.  G.  Khachian.  Polynomial  Algorithms  in  Linear  Programming.  Zhumal  Vychislitelnoi 
Matematiki  i  Matematicheskoi  Fiziki,  20:53-72,  1980. 

[7]  P.  Klein,  S.  A.  Plotkin,  C.  Stein,  and  E.  Tardos.  Faster  Approximation  Algorithms 
for  the  Unit  Capacity  Concurrent  Flow  Problem  with  Applications  to  Routings  and 
Finding  Sparse  Cuts.  Technical  Report  961,  School  of  ORIE,  Cornell  University,  1991. 

[8]  T.  Leighton,  F.  Makedon,  S.  Plotkin,  C.  Stein,  E.  Tardos,  and  S.  Tragoudas.  Fast 
Approximation  Algorithms  for  Multicommodity  Flow  Problems.  In  Proc.  23st  Annual 
ACM  Symposium  on  Theory  of  Computing,  1991. 

[9]  J.  B.  Orlin.  A  Faster  Strongly  Polynomial  Minimum  Cost  Flow  Algorithm.  In  Proc. 
20th  Annual  ACM  Symposium  on  Theory  of  Computing,  pages  377-387, 1988. 

[10]  S.  A.  Plotkin,  D.  Shmoys,  and  E.  Tardos.  Fast  Approximation  Algorithms  for  Fractional 
Packing  and  Covering.  In  Proc.  32nd  IEEE  Annual  Symposium  on  Foundations  of 
Computer  Science,  1991.  To  appear. 

[11]  F.  Shahrokhi  and  D.  Matula.  The  Maximum  Concurrent  Flow  Problem.  J.  Assoc. 
Comput.  Mach.,  37:318-334,  1990. 

[12]  P.  M.  Vaidya.  Speeding  up  Linear  Programming  Using  Fast  Matrix  Multiplication.  In 
Proc.  30th  IEEE  Annual  Symposium  on  Foundations  of  Computer  Science,  1989. 


11 


